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ABSTRACT 

We investigate the global pulsation properties of DBV white dwarf models that include both 
the double-layered envelope structure expected from time-dependent diffusion calculations, 
as well as a non-uniform C/O core expected from prior nuclear burning. We compare these 
models to otherwise identical models containing a pure C core to determine whether the ad- 
dition of core structure leads to any significant improvement. Our double-layered envelope 
model fit to GD 358 that includes an adjustable C/O core is significantly better than our pure 
C core fit (7<7 improvement). We find a comparable improvement from fits to a second DBV 
star, CBS 1 14, though the values of the derived parameters may be more difficult to reconcile 
with stellar evolution theory. We find that our models are sy stematically cooler by 1,900 K 
relative to the similar models of iFontaine & Brassard! (120021) . Although a fit to their model 
reproduces the mass and envelope structure almost exactly, we are unable to reproduce the 
absolute quality of their fit to GD 358. Differences between the constitutive physics employed 
by the two models may account for both the temperature offset and the period residuals. 
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1 INTRODUCTION 

The vast majority of what we can learn about stars is derived from 
observations of the thin outermost skins of their surface layers. For 
the most part, we are left to infer the properties of the interior based 
upon our current best understanding of the constitutive physics. 
There are a number of exceptions to this general rule, including 
observations of supernovae explosions and interacting binary stars, 
where the deeper layers of the interior are either ejected from the 
system or are gradually exposed as mass is transferred from one 
star to the other, respectively. However, both of these processes are 
disruptive to the stellar interior, making inferences about the origi- 
nal structure ambiguous at some level. 

Pulsating stars represent the best opportunity for probing stel- 
lar interiors while they are still intact. The most dramatic example 
is the Sun, where observations of light and radial velocity variations 
across the visible surface have led to the identification of thousands 
of unique pulsation modes, each sampling the solar interior in a 
slightly different and complementary way. These observations have 
led to such precise constraints on the standard solar model that the 
inverted radial profile of the sound speed, for example, now agrees 
to better than a few parts per thousand over 90 per cent of the solar 
radius IChristensen- Dalsgaarj2003) . 

If we were to move the Sun to the distance of even the nearest 
star, most of the pulsation modes that we now know to be present 
would be rendered undetectable. We would lose nearly all of our 
spatial resolution across the disc of the star, and only those modes 



stars: individual (GD 358, CBS 1 14) - stars: interiors - stars: 



of the lowest spherical degree (t < ~3) would produce significant 
variations in the total integrated light or the spectral line profiles. 
This would reduce the number of detectable modes from thousands 
to merely dozens, leading to a corresponding reduction in the abil- 
ity of the observation s to constrain the internal structure (e.g., see 
Kiel dsenetai] ri999i. Even so, such data would still allow us to 
determine the global properties of the star and to probe the gross 
internal composition and structure, providing valuable independent 
tests of stellar evolution theory. 

Pulsati ons in white dw arfs were first discovered in the cooler 
DA stars bvlLandollI fl968h . and la ter in the hotter PG 1159 stars 
iMcGraw et alll979h and DB stars IWinget et alll982t) . These dif- 
fer from the solar oscillations in that they are non-radial g-modes 
instead of p-modes; the restoring force is gravity rather than pres- 
sure, and so the most important physical quantity is the buoyancy 
frequency, not the acoustic frequency. White dwarf stars are the 
end-points of stellar evolution for all stars with masses below what 
is necessary to produce elements much heavier than carbon and 
oxygen. Their interior structure contains a record of the physical 
processes that operate during the later stages in the lives of most 
stars, so there is a potential wealth of information encoded in their 
pulsation frequencies. 

Much of the observational data for white dwarf asteroseis- 
mology has come from an interna tional c ollaboration known as 
the Whole Earth Telescope (WET: iNather et aljll990h . This net- 
work of small telescopes situated around the globe obtains the 
long time-series of photometric measurements that are often nec- 
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Figure 1. The pulsation periods observed in GD 358 (solid points) plotted 
against their deviations from the mean period spacing (dP), along with two 
physically distinct model fits (open points). The fit of Metcalfe 1 2003) has 
extra s tructure in the core (top panel), while the fit of|Fontaine & B rassard 
1 2002) has extra structure in the envelope (bottom panel). 



Figure 2. Theoretical white dwarf internal oxygen profiles from the calcula- 
tions of lSalaris eta?]<1997t dotted) and from lAlthaus et aljfeool dashed), 
along with the simplified generic profiles that we adopt for the models in 
this paper (solid) in which both the central oxygen mass fraction (Xq, ver- 
tical axis) and the fractional mass of the initial break from a uniform C/O 
mixture (q, horizontal axis) are adjustable parameters. 



essary to unambiguously identify and resolve dozens of unique 
pulsation modes. One of the most notable successes of this col- 
laboration was an observing campaign in 1990 on the brightest 
known DB variable, GD 358 JWinget et ai]ll994t) . Among many 
interesting results, these observations established the simultaneous 
presence of 11 low-degree (1=1, m=0) modes of consecutive ra- 
dial overtone (fc= 8-18) with periods in the range 400-800 seconds 
( Brad ley & Wingetll994l) . 

The theoretical interpretation of these data has grown gradu- 
ally more sophisticated as our computational capabilities have ex- 
panded over the decade following the observations. The most recent 
analysis produced a model that leads to a root-mean-square differ- 
ence between the observed and calculated periods of only ~1 sec- 
ond ^Metcalfe 2003). However, a comparabl e match to the same ob - 
servations was recently published bv lFontaine & Brassard 1 2002), 
who used a completely independent model with an inter nal struc- 
ture that was physicall y distinct from that a ssumed bv iMetcalfd 
(see Fig.0. Essentiallv. lFontaine & Brassarof s model contains ad- 
ditional structure in the envelope, while Metcalfe's model has extra 
structure in the core. 

While this duality in model-fits is now understood to be due to 
an inh erent symmetry in the way the pulsations sample the interior 
I Mont gomery. Metcalfe. & W inget 2003), there is good reason to 
believe that the physical basis of each model is sound, but that nei- 
ther represents a complete description of the true interior structure. 
In this paper, we attempt to bridge the gap between these two mod- 
els by including the essential elements of each into a 'hybrid' model 
that contains both the double-layered envelope structure expected 
from time-dependent diffusion calculations, and an adjustable car- 
bon/oxygen (C/O) core. In ij2]we outline the physical basis of our 
model parametrizations, and in J3]we present the results of our 
model-to-model and model-to-data comparisons. We conclude in 
^4]with a brief discussion of the promise of asteroseismology. 



2 THEORETICAL MODELS 
2.1 Adjustable C/O Cores 

The generic shape of a theoretical white dwarf internal oxygen 
profile is set by the nuclear and mixing processes that occur dur- 
ing its formation in the core of a red giant. The detailed shape is 
less certain, since it depends on the specific physical and numerical 
treatments utilized by the adopted model. Many models agree that 
the inner ~0.5 in fractional mass should contain an approximately 
uniform C/O mixture JSalaris et alll997tlAl7haus et all2002T) . The 
precise extent of this region is determined by the maximum size of 
the central convective core (and by the amount of convective over- 
shooting, if it occurs) during helium burning. The C/O ratio in this 
uniform region is primarily set by the rate of the 12 C(a,7) 16 re- 
action. Outside of the uniform C/O core the exact manner in which 
the oxygen mass fraction drops to zero is uncertain, since it de- 
pends on the adopted recipe for semiconvection during core helium 
burning and on the physical conditions in the models during helium 
shell burning. However, the location and slope of the initial break 
from a constant C/O ratio seems to be the most important feature 
from an asteroseismological standpoint. 

Our parametrization for the C/O core is iden- 
tical to that us ed in the earlier mod el fitti ng of 
I Metcalfd 120031). IMetcalfe. Salads. & Winged §002), and 
iMetcalfeWmgeT^Toiarb^ fix the oxygen mass 

fraction to its central value (Xo) out to some fractional mass (q) 
where it then decreases linearly in mass to zero oxygen at the 
0.95 fractional mass point (see Fig.|2j- This form for the profile 
was chosen to facilitate comp arison with the earlier work of 
iBradlev. Winget. & Wood! 19931) 

Despite this simplification, IMetcalfe! J2003) obtained fits to 
two different pulsating DB white dwarfs that independently led to 
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Figure 3. Internal He pr ofiles from the cal culations of Dehn er & Kawalei 
1 1995, dotted) and from lFontaine & Brassard] <200l dashed), along with a 
profile from the parametrizafion that we adopt for the models in this paper 
(solid). The locations of the two chemical transitions (points) are adjustable 
parameters, but the shapes are fixed. 

implied rates for the 12 C(a, 7) 16 reaction that were consistent both 
with each other and with recent extrapolations from high energy 
laboratory measurements. While these fits are the best evidence to 
date that our models really are sensitive to the core composition 
and structure sampled by the white dwarf pulsations, they assumed 
a single-layered structure for the transition between the C/O core 
and the pure He envelope. If there is an evolutionary connection 
between the DB stars and their hotter cousins, the PG 1159 stars, 
then this single-layered structure may not be an adequate represen- 
tation of the actual white dwarf envelopes. 

2.2 Double-Layered Envelopes 

iDehner & Kawaletl Il995l) were the first to demonstrate the plau- 
sibility of an evolutionary connection between PG 1159 and DB 
stars. They c reated an initial m odel based on the pulsation data 
for PG 1 159 Iwinget et alll99ll) . which included a thick envelope 
[logfMcv/M*) ~ -2.5] containing a uniform mixture of 30% He, 
35% C, and 35% O above the C/O core. During the evolution of this 
model to lower temperatures, they included time-dependent diffu- 
sive processes so they could monitor the chemical and gravitational 
separation of the elements - essentially He 'floating' to the surface. 
By the time the model had reached the temperature range of the DB 
variables, the transition between the pure He surface and the still- 
uniform He/C/O mantle was located near log(MH e /M») ~ -5.5, 
leading to a structure that looked very similar to the outer regions 
of the fit of Bradley & Winget ( 1994). The pulsation periods of this 
model were a good match to the WET observations of GD 358, 
even though no formal attempt w as made to fit t hem. 

More recently, Fontaine & Brassard (2002) calculated a large 
grid of similar models that included this double-layered envelope 
structure in an attempt to match the pulsation periods of GD 358 in 
detail. Using models that contained a pure C core, they were able to 



find a matc h with root-mea n-square period residuals nearly as low 
as the fit of Metcalfe 1 2003). They also attempted similar fits using 
models with cores of pure O and a uniform C/O mixture, resulting 
in models that were only 'marginally worse'. They speculated that 
the 1 1 periods observed in GD 358 may not contain enough infor- 
mation to measure the core and envelope structure simultaneously. 

To investigate this possibility more thoroughly, we replaced 
the single-layered envelopes in our models with a p arametriza- 
tionof the double-layered structure used by iFontaine & Brassard! 
(2002). This allows us to specify the locations in the enve- 
lope where the He mass fraction reaches 0.15 (the middle of 
the transition between the C/O core and the He/C mantle) and 
0.65 (the middle of the transition between the mantle and the 
pure He surface) 1 . The He profiles from several different mod- 
els are shown in Fig. [3] Although there are slight differences be- 
tween the profile from our parameterization compared to that of 
Fonta ine" & Brassardl 120021) . we have verified that modifications 
to the detailed shape of the profile comparable to these differ- 
ences lead to only marginal improvements in the final fits. Note 
that the two layers are theoretically expected to merge into one at 
low values of 7rff l Gautschv & Althaus 1 l2002l) . leading to struc- 
tures that have effecti vely been explored by the sin gle-layered 
glo bal model-fitting of iMe tcalfe. Winget. & CharbonneaJ (2001) 
and lHandler. Metcalfe. & Wood! 120021) for GD 358 and CBS 114 
respectively. 

3 NEW MODEL FITTING 

To get a sense of the relative improvement that might be possible 
by including an adjustable C/O core in our models, we performed 
two sets of fits to the puls ation data. In both cases, w e used a ge- 
netic algorithm (see Metcalfe & Charbonneau 2003) to minimize 
the root-mean-square residuals between the observed and calcu- 
lated periods (op) for models with effective temperatures (Jkf) be- 
tween 20,000 and 30,000 K, and total stellar masses (Af*) between 
0.45 and 0.95 Mq, We allowed the location of the inner He tran- 
sition (between the core and the mixed He/C mantle) to assume 
values of log(M env /M») between -2.0 and -4.0, and for the outer 
He transition (between the mixed He/C mantle and the pure He 
surface) the value of log(MH e /M») was in the range -5.0 to -7.0. 
For the first set of fits we fixed the core composition to be pure C. 
For the second set we included an adjustable C/O core, allowing 
the central oxygen mass fraction (Xo) to vary between 0.0 and 1.0, 
with the break from a uniform C/O mixture beginning at a frac- 
tional mass (q) between 0.10 and 0.85. 

We performed these fits on three sets of pulsation data. To ex- 
plore the systematic differences between our own models and those 
of Fon taine & Br assard ( 200^), we fit the 1 1 model pulsation peri- 
ods of their fit to GD 358, which had a pure C core. For our ad- 
justable C/O fit to these periods, we used a single-layered envelope 
model as part of an investigatio n of the symmetry bet ween these 
two types of models (Montgomery, Metcalfe, & Winget 2003). The 
second set of data came from the 1990 WET campaign on GD 358 
( Winge^yriJ^}94|), which is the same set of 11 periods fit by 
Fontaine & Brassard (2002). To investigate whether the double- 
layered envelope models could produce fits to a second DBV star 
that were physically consistent with time-dependent diffusion the- 
ory, our third set of data came from the slightly hotter white 

' IFonta ine & Bras sard define the location of the outer transition to be 
where the He mass fraction reaches 0.50, which is systematically deeper 
by ~0.3 dex relative to our definition. 
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Table 1. Double-layered envelope models, with and without an adjustable C/O core, fit 
to the model periods of Fontaine & Brassard 1 2002 ) and to the periods observed in the 
pulsating DB white dwarf stars GD 358 and CBS 1 14. 



F&B Model GD 358 CBS 114 

Parameter pure C C/O" pure C C/O - pure C C/O 



7 e ff (K) 22,900 22,900 21,500 20,800 24,500 25,800 

M* (Mq) 0.625 0.660 0.665 0.670 0.650 0.615 

log(M e „v/M„)... -2.94 -2.00 -3.10 -2.68 -3.96 -3.78 

log(M Hc /M»)... -6.10* ••• -6.22 -5.22 -5.72 -5.74 

X ■•• 0.67 ■■■ 0.75 ■■• 0.71 

q ■•• 0.48 ■•• 0.49 ■■• 0.56 

crp(s) 1.63 1.26 2.17 1.52 0.82 0.27 



a Single-layered envelope fit from Montgomery. Metcalfe. & Wingej] 120031) 
^ By Fontaine & Brassard's definition, log(MH c /M») = -5.77 



dwarf CBS 114, whic h exhibits 7 independent pulsation periods 
(Handler, Metcalfe. & Wood 2002). In all cases we restricted our 
models to calculate the £=l, m=0 modes, as suggested by the ob- 
servations. The results of our two fits to each of these three data 
sets are shown in TableQ 

3.1 Carbon Core Models 

From our pure C core fit to the model periods of 
iFontaine & Brassarcl 12002L their table 1) we retrieved the 
same mass and virtually identical locations for the two He transi- 
tions in the envelope, but our temperature was systematically low 
by 1,900 K relative to their value (24,800 K). Because the internal 
structure of the two models is very similar, this temperature offset 
is most likely due to differences between the constitutive physics 
(equations of state, opacities, etc.) as well as the convective 
prescription employed by each. In a study of the influence of 
various updates to the cons titutive physics on the pul sational 
properties of their models, Fontaine & Brassard 1 1994) found 
that the OPAL radiative opacities faglesias^^RogersI 1 19931) for 
pure He and C were syst ematically lower than the older LAO 
data iHuebner et alj Il97?1) . and that models using the higher 
opacity data could mimic hotter models with otherwise identical 
parameters. This is consistent with the result s of our fit, since our 
models use the more opaque LAO data, while Fontaine & Brassard 
use th e more recent OPAL data (Fontaine, Brassard. & Bergeron 
1200 ll) . However, our preliminary attempts to use the OPAL 
data in our own code suggest that this may only account for 
roughly half of the temperature difference between the models. 
Other possible contributions include differences between our 
equati ons of state - both of our models use tables derived from 
Lamb & Van Hornl 1 19751) and IFontaine Graboske & Van Horn! 
Il977l) for p arts of th e interior and the envelope respectively, but 
Font aine" & Brassards models include additional unpublished 
modifications iFontaine, Brassard, & Bergeron 2001). More work 
will be required to document these differences and to evaluate their 
influence on the effective temperature scales of our models. 

Our fit to the observed pulsation periods of GD 358 
has a higher mass and a lower temperature than the fit of 
IFontaine & Brassarcl l2002h . but the locations of the two He tran- 
sitions are both within their la uncer tainties. The most st riking 
difference between our fits is that while IFontaine & Brassarcl were 
able to match the periods observed in GD 358 with residuals of 
only <rp = 1 .30 s, our optimal fit has at = 2. 17 s (a difference nearly 



30 times larger than the observational uncertainty, a a u ~ 0.03 s). 
Again, differences between the constitutive physics employed by 
the two models may account for at least part of this difference in the 
absolute quality of the two fits - using the OPAL radiative opacities 
with our models led to a significant decreas e in the residuals (e.g ., 
from 1.63 to 1.45 s in the case of our fit to IFontaine & Brassarot s 
model). Despite the unresolved differences between the constitu- 
tive physics, we can still use the relative quality of our own fits to 
determine whether or not the addition of extra free parameters in a 
given model lead to significant improvements. 

Our double-layered envelope fit to CBS 1 14 is only marginally 
(2<7o bs) better than the pure C core sing le-layered envelope model 
fit of Handler, Metcalfe, & Wood (2002), compared to what is ex- 
pected from the addition of an extra parameter. Our values for the 
mass, the temperature, and the (inner) He transition are indistin- 
guishable from the single-layered fit. Although the derived temper- 
ature for CBS 1 14 is higher than for GD 358, as expected from the 
spectroscopic measurements of iBeauchamp et al] 1 19991) . the pure 
He surface layer for this fit is thicker than for GD 358 - just the 
opposite of what time-dependent diffusion calculations would lead 
us to expect for models with similar masses. 

3.2 C/O Core Models 

When we used an adjustable C /O core and a single- l ayered enve- 
lope to fit the model periods of IFontaine & Brassarcl 120021) . it led 
to a significant improvement to our match even though the interior 
structure of our model was dramatically different from the source 
model. Our derived value of q corresponds exactly to the true loca- 
tion of the outer He tra nsition refl ected thr ough the co re/envelope 
symmetry described by M ontgomery. Metcalfe. & Wingetl 120031) . 
Essentially our fit confirms empirically that such a symmetry ex- 
ists, and that it is possible to fit real structure in the envelopes by 
assuming structure at the corresponding symmetric location in the 
core. 

Our fit to GD 358 leads to inferred locations for all chem- 
ical transition zones that are distinct from each other with re- 
spect to the core/envelope symmetry. Statistically, the addition of 
an adjustable C/O core to the double-layered envelope fit leads 
us to expect the residuals to decrease to 1.74 s just from the ad- 
dition of the two extra parameters. In fact the residuals of our 
C/O fit are reduced to 1.52 s, a 7o- b s improvement beyond what 
is expected. In this model, the C/O core is effectively fitting the 
same mode trapping structure attributed to the outer He transition 
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in iFontaine & Brassard ' s fit. Our outer He transition, centered at 
log(MHe/M„) = -5.22, is distinct since no structure is expected at 
the corresponding symmetric locatio n in the core, which is near 
Mr/M* ~ 0.6 (see Montgomery. Metcalfe, & Winaet 2003, their 
fig. 3). 

The improvement in the fit to CBS 114 is even more sub- 
stantial. While the derived locations of the two He transitions do 
not change by much between the pure C and the C/O fits, the 
addition of structure in the core leads to an improvement in the 
fit of nearly 12<T b s beyond what is expected from the extra pa- 
rameters. As with GD 358, the locations of the core and outer 
envelope structure are distinct from one another, from the stand- 
point of the symmetry inherent in the models. Interestingly, the 
conflict between the higher temperature and the thickness of the 
surface He layer disappears for the C/O fits. However, both fits 
to CBS 114 lead to an unusually small total envelope mass rel- 
ative to the range expected [log(M e „v / M*) ~ -2.0 to -3.0] from 
simulations of carbon dredge-up in DQ stars (see IPelletier et all 
Il986t IFontaine & Brassard KOfO. In addition, the masses of the 
C/O fits are in c onflict with the spectroscopic measurements of 
iBeauchampetalH 19991) . 

4 DISCUSSION 

Our global exploration of models containing double-layered enve- 
lope structure make it clear that the addition of an adjustable C/O 
core leads to significantly better fits to the observations, relative to 
pure C core models. This is reassuring, since there are sound phys- 
ical reasons that lead us to expect composition transition zones in 
both the cores and the envelopes of real white dwarf stars. 

The absol ute quality of our fits to GD 358 are generally wo rse 
than the fits of IFontaine & Brassard l2002t) or lMetcalfd 120031) . If 
white dwarf envelopes really do contain a double-layered structure, 
this implies that we will need to know the detailed shapes of the 
composition transition zones in both the core and the envelope be- 
fore better fits can be found. 

The fits to CBS 114 support our general conclusion that the 
models are sensitive to structure in both the core and the envelope, 
but only the C/O fit is consistent with the expectations from time- 
dependent diffusion theory (thinner surface He layers for hotter 
white dwarfs of comparable mass). However, the relatively small 
derived values for the stellar mass and the total envelope mass re- 
main a concern, implying that double-layered models may be less 
able to explain the observations of both GD 358 and CBS 114 in 
a physically self-consistent manner (cf. Metcalfe 2003). Additional 
calculations of the mass-dependence of time-dependent diffusion 
profiles in white dwarf envelopes and higher signal-to-noise obser- 
vations of CBS 114 to yield a larger number of periods for fitting 
would both clearly be useful. 

Comparison of our model to that of Fontain e* & Brassard" 
revealed a systematic difference of 1,900 K between our ef- 
fective temperature scales. This may help us to understand the 
relatively low effective temperatures derived in previous stud- 
ies using our models "M etcalfe. Winget, & Charbo nnead 1200 it 
Metcalfe, Salaris, & Wingetl l2002t lMetcalfdl2003l) . The source of 
this temperature offset appears to be related to our use of the larger 
LAO radiative opacities, and possibly due to subtle differences be- 
tween our convective prescriptions and the equations of state that 
we use for both the core and the envelope. 

Our understanding of stellar interiors through white dwarf as- 
teroseismology is evolving rapidly. It is now clear that the pul- 
sations really do sample the star globally, and we must be care- 



ful to avoid the potential ambiguities caused by the intrinsic 
core/envelope symmetry. Our challenge is to find a physically self- 
consistent description of both the core and the envelope that will 
allow us to match the pulsation periods within the observational un- 
certainties. With persistence, and with an open-source development 
philosophy for our models ( Metc alfe et"aTl2 002). the future holds 
great promise for unveiling the detailed composition and structure 
of white dwarf interiors. 
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